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We apply the postquasistatic approximation, an iterative method for the evolution of self- 
gravitating spheres of matter, to study the evolution of anisotropic non-adiabatic radiating and dis- 
sipative distributions in General Relativity. Dissipation is described by viscosity and free-streaming 
radiation, assuming an equation of state to model anisotropy induced by the shear viscosity. We 
match the interior solution, in non-comoving coordinates, with the Vaidya exterior solution. Two 
simple models are presented, based on the Schwarzschild and Tolman VI solutions, in the non- 
adiabatic and adiabatic limit. In both cases the eventual collapse or expansion of the distribution 
is mainly controlled by the anisotropy induced by the viscosity. 

PACS numbers: 04.25.-g,04.25.D-,0.40.-b 



I. INTRODUCTION 

In order to study astrophysical fluid dynamics one can 
get complicated models incorporating realistic transport 
mechanisms and equations of state. The simplest case 
with mass, spherical symmetry, despite its simplicity, 
still remains an interesting problem in numerical rela- 
tivity, specially when including dissipation. Dissipation 
due to the emission of massless particles (photons and/or 
neutrinos) is a characteristic process in the evolution of 
massive stars. It seems that the only plausible mecha- 
nism to carry away the bulk of the binding energy of the 
collapsing star, leading to a black hole or neutron star, 
is neutrino emission Viscosity may be important in 
the neutrino trapping during gravitational collapse 0-01 , 
which is expected to occur when the central density is of 
the order 10^^-lO^^g cm^'^. Although the mean free path 
of the neutrinos is much greater than for others particles, 
the radiative Reynolds number of the trapped neutrinos 
is nevertheless small at high density Q , rendering viscous 
the core fluid 0, [H. 

Numerical Relativity is expected to keep its power to 
solve problems and generate new, interesting physics, 
when dissipative distributions of matter are considered. 
In fact, numerical methods in General Relativity have 
been proven to be extremely valuable for the investi- 
gation of strong field scenarios (see and references 
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therein). For instance, these methods and frameworks 
have (i) revealed unexpected phenomena [1] , (ii) enabled 
the simulation of binary black holes (neutron stars) [olIToj 
and (iii) allowed the development of relativistic hydrody- 
namic solvers [ll|, among other achievements. Currently, 
the main limitation for Numerical Relativity is the com- 
putational demand for 3D evolution . The addition of 
a test-bed for studying dissipation mechanisms and other 
transport processes in order to later incorporate them 
into a more sophisticated numerical framework (ADM or 
characteristic) is a necessity. 

In this paper, we study a selfgravitating spherical dis- 
tribution of matter containing a dissipative fluid. We 
follow the method proposed in [l^, which introduces a 
set of conveniently defined "effective" variables (effective 
pressure and energy density), where their radial depen- 
dence is chosen on heuristic grounds. In essence this is 
equivalent to going one step further from the quasistatic 
regime, and the method has been named the postqua- 
sistatic approximation (PQSA) after [3 • The essence 
of the PQSA was first proposed in [l^ using radiative 
Bondi coordinates and it has been extensively used by 
Herrera and collaborators (T6l - [23| . By quasistatic ap- 
proximation we mean that the effective variables coin- 
cide with the corresponding physical variables (pressure 
and energy density). However, in Bondi coordinates the 
notion of quasistatic approximation is not evident: the 
system goes directly from static to postquasistatic evolu- 
tion. In an adiabatic and slow evolution we can catch-up 
that phase, clearly seen in non-comoving coordinates. 
This can be achieved using Schwarzschild coordinates 
Here we study radiating viscous fluid spheres in 



2 



the streaming out limit with the PQSA approach which 
allow us to departure from equilibrium in non-comoving 
coordinates. These systems have been studied using the 
method described in Ts] for the radiative shear viscos- 
ity problem and its effect on the relativistic gravitational 
collapse [IS-dll- We do not consider temperature pro- 
files to determine which processes can take place during 
the collapse. For that purpose, transport equations in 
the relaxation time approximation have been proposed 
to avoid pathological behaviors (see for instance [22| and 
references therein). These issues will be considered in 
a future investigation. In order to develop a numerical 
solver which incorporates in a realistic way dissipation 
following the Miiller-Israel-Stewart theory [29l - l32l | it is 
first necessary to know, to zeroth level of approximation, 
viscosity profiles like the ones presented in this investiga- 
tion. The physical consequences of considering dissipa- 
tion by means of an appropriate causal procedure have 
been stated analytically in several pap ers by Herrera and 
collaborators (see for example |33l-l36{). 

To the best of our knowledge, no author has under- 
taken in practice the dissipative matter problem in nu- 
merical relativity. Our purpose here is to show how vis- 
cosity processes can be considered as anisotropy and how 
the PQSA works in this context. Our results partially 
confirm previous investigations psl [26| . The novelty here 
is in the use of the PQSA to study dissipative scenarios. 
The results indicate that an observer using radiation co- 
ordinates does not " see" some details when shear viscos- 
ity is considered. The final goal is to eventually study the 
same problem using the Miiller-Israel-Stewart theory for 
dissipative system, which is highly nontrivial in spherical 
symmetry. 

In standard numerical relativity, in order to deal with 
matter both in ADM 3-1-1 [37] and in the characteristic 
formulations [11], Bondian observers have been used im- 
plicitly. This has been noted recently, and the method 
has been proposed as a test bed in numerical relativity 
[38| . The systematic use of local Minkowskian and co- 
moving observers in the PQSA, named Bondians, was 
used to reveal a central equation of state in adiabatic 
scenarios [s^, and to couple matter with radiation [40| . 
Since Bondian observers are a fundamental part of the 
PQSA and all its applications in the characteristic for- 
mulation we are currently trying to transfer all the expe- 
rience gained using this approach to include more real- 
istic effects in the dynamics of the fluid using the ADM 
3+1 formulation, the most popular method in numerical 
relativity. Besides introducing a more realistic time scale 
in the problem with matter the intention is to promote 
the PQSA (and any of its applications) as a test bed in 
the ADM 3-1-1 and characteristic approaches. 

The plan for this paper is as follows: In Section ^|TT] 
we present the field equations and matching conditions 
at the surface of the distribution. We explain the PQSA 
and write a set of surface equations, in Section ijllll In 
Section jjlVl we illustrate the method presenting four sim- 
ple models based on the Schwarzschild and Tolman VI in- 



terior solutions. Finally, we discuss the results in Section 



II. FIELD EQUATIONS FOR BONDIAN 
FRAMES AND MATCHING 



To write the Einstein field equations we use the line 
element in Schwarzschild~like coordinates 

ds^ = e^dt^ - e^dr^ - {d9^ + sin ^0d,p^) , (1) 

where v — v{t^r) and A = \{t,r), with (t,r,0,(f>) = 
(0,1,2,3). 

In order to get physical input we introduce the 
Minkowski coordinates (t, x, y, z) by [41] 

dr = e'^/^di, dx = e^/'^dr, dy = rd0, dz = r sin 6*^0, (2) 

In these expressions v and A are constants, because they 
have only local values. 

Next we assume that, for an observer moving relative 
to these coordinates with velocity uj in the radial (x) 
direction, the space contains 

• a viscous fluid of density p, pressure p, effective bulk 
pressure and effective shear pressure p^, and 

• unpolarized radiation of energy density e. 

For this moving observer, the covariant energy tensor 
in Minkowski coordinates is thus 



p + e -e 

-e P + f — PC — '^Pr) 

P-Pc+P,, 

P-Pc+Pr, 

(3) 

Note that from ^ the velocity of matter in 
Schwarzschild coordinates is 



dr 
'di 



(4) 



Making a Lorentz boost and defining p = p — p^, Pr = 
p — 2pri, pt = P + Pjj and e = e(l -I- w)/(l — uj) we write 
the field equations in relativistic units (G = c = 1) as 
follows: 
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where the comma (,) represents partial differentiation 
with respect to the indicated coordinate and the effec- 
tive variables p, 5, known as conservation variables as 
well, and p, the flux variable, 
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S = (p + Pr 



and 



p = 



Pr + pUJ 
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Equations ([S])-® are formally the same as for an 
anisotropic fluid in the streaming out approximation [2^ . 

At this point, for the sake of completeness, we write the 
effective viscous pressures in terms of the bulk viscosity 

the volume expansion Q, the shear viscosity rj and the 
scalar shear a 



and R is a null coordinate {gRR = 0). The relation- 
ship at the surface between the coordinates {t,r,9,(j)) and 
{u,R,9,(f>) is 



u = t - r — 2M In 



2M 



-l],R = r. 



(17) 



The exterior and interior solutions are separated by 
the surface r — a{t). In order to match both regions 
on this surface we use the Darmois junction conditions. 
Demanding the continuity of the first fundamental form, 
we obtain 
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From now on the subscript a indicates that the quantity 
is evaluated at the surface. Now, instead of writing the 
junction conditions as usual, we demand the continuity 
of the first fundamental form and the continuity of the 
independent components of the energy-momentum flow 
(4^ . This last condition guarantees the absence of sin- 
gular behaviors on the surface. It is easy to check that 



PC = C©, 



(12) 



Pa = PCa + '^PVa 



(20) 



Pn - ^W, 



where 



e = 



(1 -a;2)i/2 
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and 
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We have four field equations for five physical variables 
(p, Pr, e, id and pt) and two geometrical variables (v and 
A). Obviously, we require additional assumptions to han- 
dle the problem consistently. However, we discuss first 
the matching with an exterior solution and the surface 
equations that govern the dynamics. 

We describe the exterior space-time by the Vaidya 
metric 
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2M{u) 
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du^+2dudR~R^ 
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(16) 

where w is a time-like coordinate so that u = constant 
represents, asymptotically, null cones open to the future 



which expresses the discontinuity of the radial pressure 
in the presence of viscous processes. 

Before proceeding with the description of the method 
it is convenient to rewrite some equations and introduce 
one equation of state. 

Defining the mass function as 



e"^ = 1 - 2m/r, 

and substituting ([2T|l into ([5]) and | 
some arrangements, 



(21) 

we obtain, after 



dm 
~dt 



^ + e(l-t^)(l-2m/r)i/2g^/2 



This equation, known as the momentum constraint in the 
ADM 3+1 formulation, expresses the power across any 
moving spherical shell. 

Equation ([7|) can be written as T^.^^ = or equivalently, 
after a lengthly calculation 

~ , ip + p){47Tr^p + m) 2 
r[r — Zm) r 



47rr(r — 2m) 



m,tt 



" — 2m 



(23) 



This last equation corresponds to a generalization of 
the hydrostatic support equation, that is, the Tolman- 
Oppenheimer-Volkoff (TOV) equation. It can be shown 



4 



that equation ((23)) is equivalent to the equation of mo- 
tion for the fluid in conservative form in the standard 
ADM 3+1 formulation Equation leads to the 
third equation at the surface (see next section); up to this 
point is completely general within spherical symmetry. 

To close this section we have to mention that we 
assume the following equation of state for non- 
adiabatic modeling [43] 



Pt -Pr 



C{p + p){4:'Kry + m) 
(r — 2m) 



(24) 



where C is a constant. 
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FIG. 1: Evolution of the radius A{t) for the Schwarzscliild- 
like model I. The initial conditions are ^4(0) — 5.0, -F(O) = 0.6, 
n{0) = -0.1. 
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FIG. 2: Evolution of the energy density p (multiplied by 10'') 
for the Schwarzschild-like model I. The initial conditions are 
A{0) = 5.0, F{0) = 0.6, 0.(0) = -0.1, with h = 0.99. 
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III. THE POSTQUASISTATIC 
APPROXIMATION AND THE SURFACE 
EQUATIONS 

Feeding back © and ((TTl) and using ([H]) into (O and 
, these two field equations may be formally integrated 
to obtain 

m = / Anr'^p dr (25) 

which is the Hamiltonian constraint in the ADM 3+1 
formulation and 

r 2(47rr3p + m) , 
J„ r(r - 2to) 

the polar slicing condition, from where it is obvious that 
for a given radial dependence of the effective variables, 
the radial dependence of the metric functions becomes 
completely determined. 
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FIG. 3: Evolution of the radial pressure pr (multiplied by 10"^) 
for the Schwarzschild-like model I. The initial conditions are 
A{Q) = 5.0, F{Q) = 0.6, 0(0) = -0.1, with h = 0.99. 



As defined in [ij] the postquasistatic regime is a sys- 
tem out of equilibrium (or quasiequilibrium; see |34|) 
but whose effective variables share the same radial de- 
pendence as the corresponding physical variables in the 
state of equilibrium (or quasiequilibrium). Alternatively, 
we can say that the system in the postquasistatic regime 
is characterized by metric functions of the static (qua- 
sistatic) regime. The rationale behind this definition is 
not difficult to catch: we look for a regime which, al- 
though out of equilibrium, it is the closest to quasistatic 
evolution. 
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FIG. 4: Evolution of the tangential pressure pt (multiplied by 
10^) for the Schwarzschild-like model I. The initial conditions 
are A{0) = 5.0, F{0) = 0.6, 0(0) = -0.1, with h = 0.99. 



FIG. 6: Evolution of the radiation flux e (multiplied by 10*) 
for the Schwarzschild-like model I. The initial conditions are 
A{0) = 5.0, F{0) = 0.6, 0.(0) = -0.1, with h = 0.99. 
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FIG. 5: Evolution of the local velocity lj for the 
Schwarzschild-like model I. The initial conditions are A{0) — 
5.0, F(0) = 0.6, fl(0) = -0.1, with h = 0.99. 



FIG. 7: Evolution of the expansion O for the Schwarzschild- 
like model I. The initial conditions are ^4(0) = 5.0, F{0) = 0.6, 
n{0) = -0.1, and h = 0.99. 



A. The PQSA protocol 

We outline here the PQSA approach: 

1. Take an interior solution to Einstein's field equa- 
tions, representing a fluid distribution of matter in 
equilibrium, with static solutions 

Pst. = p(r), Pst. =p{r). 

2. Assume that the r dependence of p and p is the 
same as that of p^t. and Pst., respectively. 



3. Using equations (^51 and (pS)) . with the r depen- 
dence of p and p, one gets m and v up to some 
functions of t, which will be specified below. 

4. For these functions of t one has three ordinary dif- 
ferential equations (hereafter referred to as the sur- 
face equations), namely: 

(a) Equation Q evaluated at r = a; 

(b) Equation ([221) evaluated at r = a; 

(c) Equation T^.^ = evaluated at r = a. 



6 



0) 



.18 
. 16 
0.14 
. 12 
. 1 
.08 
. 06 
.04 
. 02 


-0 . 02 



r/a=0 


2 


r/a=0 


4 


r/a=0 


6 


r/a=0 


8 


r/a=l 






4 5 6 
Time 



7 8 



10 



FIG. 8; Evolution of the shear a for the Schwarzschild-hke 
model I. The initial conditions are A{0) = 5.0, F(0) = 0.6, 
n(0) = -0.1, and h = 0.99. 
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FIG. 9: Evolution of the shear viscosity rj for the 
Schwarzschild-like model I. The initial conditions are ^(0) — 
5.0, F(0) = 0.6, Q{0) = -0.1, with h = 0.99. 



7. Feeding back the result of integration in the expres- 
sions for m and v, these two functions are com- 
pletely determined. 

8. With the input from point 7, and using the field 
equations, together with the equation of state 
and/or transport equation, all physical variables 
can be found everywhere inside the matter distri- 
bution. 

As it should be clear from the above, the crucial point 
in the algorithm is the system of equations at the surface 
of the distribution. We specify it in the next section. 



B. Surface equations 



Evaluating (1^ at the surface and using the boundary 
condition ([^0]). the energy loss is given by 



-47ra^ea(l - 2ma/a){l -~ Wq). 



(27) 



Hereafter the overdot indicates d/dt and the a subscript 
indicates that quantity is evaluated at the surface r = 
ait). 

The evolution of the boundary is governed by equation 
(HJ evaluated at the surface 



a = (1 - 2ma/a)uJa- 



(28) 



Scaling the total mass ma , the radius a and the time-like 
coordinate by the initial mass ma{t = 0) = ma(0), 

A = a/ma{0), M = ma/ma{0), V"ia(0) t, 

and defining 

2M 



F= 1 



A 



n = OJa, 

the surface equations can be written as 



(29) 
(30) 
(31) 
(32) 



5. Depending on the kind of matter under consider- 
ation, the system of surface equations described 
above may be closed with the additional informa- 
tion provided by the transport equation and/or 
the equation of state for the anisotropic pressure 
and/or eventual additional information about some 
of the physical variables evaluated on the surface of 
the boundary (e.g. the luminosity). 

6. Once the system of surface equations is closed, it 
can be integrated for any initial data. 



F^^[{1- F)n + 2E] . (33) 

Equations (j32p and (j33l) are general within spherical sym- 
metry. 

We need a third surface equation to specify the dy- 
namics completely for any set of initial conditions and a 
given luminosity profile E[t). For this purpose we can 
use the field equation ([7]) or the conservation equation, 
(|23p . written in terms of the effective variables, which is 
clearly model-dependent. 
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FIG. 10: Evolution of the radius A(t) for the Schwarzschild- 
Uke model II. The initial conditions are ^(0) = 5, F{0) — 0.6, 
n{0) = -0.01. 
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FIG. 12: Evolution of the shear viscosity rj (multiplied by 10^) 
for the Schwarzschild-like model II. The initial conditions are 
A{0) = 5, F{0) = 0.6, n{0) = -0.01, and h = 0.9. 
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FIG. 11: Evolution of the local radial velocity uj for the 
Schwarzschild-like model II. The initial conditions are A(Q) = 
5, F(0) = 0.6, n{0) = -0.01, and h = 0.9. 



A. Schwarzschild-like model I: non— adiabatic 

We consider here a very simple model inspired by the 
well-known Schwarzschild interior solution 441 . We take 



fit), 



(34) 



where / is an arbitrary function of t. The expression for 
p is 



P 



P + P 



h/2 



k{t), 



(35) 



where fc is a function of t to be defined from the boundary 
condition (j20p . which now reads, in terms of the effective 
variables, as 



Pa ^ Pa^^ + ia{l + ■ 

Thus, dSSl) and §Q give 

3(1 -F) 



P = 



87ra2 
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IV. EXAMPLES - _ P [ XsF^^'^ - SV'sC 1 

^" 3 \ ^se-XsFV2 J' 

We illustrate the PQSA method with four examples -^vith 
based on the Schwarzschild and Tolman VI interior so- 
lutions. Additionally, we consider two corresponding ^ = [1 — (1 — F)(r/a)^]''^^ 
adiabatic models, that is, without free-streaming but 

anisotropic (viscous). Although greatly simplified, the where h = 1 — 2C and 
adiabatic models lead to non-trivial results, which allow 

to understand our results better. xs = 3(1^^ + 1)(1 - F) + 2E{1 + ft), 



(38) 
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ij^s = (Sfl^ + 1)(1 - F) + 2E[1 + n). 

Using ((2T|) and ((26)) it is easy to obtain expressions for 
771 and V. 



m = niair/af, 



Xs 



2(1 -F) / 



2/h 



(39) 



(40) 



In order to write down explicitcly the surface equations 
for this example, we evaluate the equation ([23]) at the 
surface, obtaining 



n 



- 8E^ 



= [8EF + lon^F - 6En + 2En^ + 
- m^F^ - en'^F + seq^ + sF^n^ + AE^n 4 

+ iEA + GFEfl - 2Fn^E - 8Fn^E) 

I (2A(F-1)) (41) 



It is interesting to note that this equation is the same 
as in the isotropic case (p,- — pt). This is a direct con- 
sequence of the chosen equation of state combined with 
incompressibility of the fluid; it is not a general result, as 
we will see for the next models. Equation (PT|). together 
with and constitute the system of differential 
equations at the surface for this model. It is necessary 
to specify one the luminosity as a function of t and the 
initial data. We choose i? to be a gaussian 



E = Eoe' 



-(t-to)Vs^ 



0.25, which 
= 1/10 of the 



with Eo = Mr/^/Sn, to = 5.0 and E = 
corresponds to a pulse radiating away Mr 
initial mass. 

We solve equations (j32|) . (|33|) . and (|4T|) using a fourth 
order Runge-Kutta method. The physical variables (p, 
p, uj, rj, e) are obtained from the field equations (l5])-(|8]) 
and the equation of state (|24|) . Note that we have to use 
equations (jl4p and (jl5p and some additional numerical 
work to determine 77. We take as initial conditions A{0) = 
5, M(0) = 1, fl{0) = -0.1, with h = 0.99. 

Figure 1 shows the evolution of the radius of the dis- 
tribution. Figures 2-6 display the physical variables (p, 
Pr, Pt, OJ, e), figures 7-8 the kinematic variables O and 
a, and figure 9 the shear viscosity rj, for different regions. 
It is evident that the emission of energy decreases the 
energy density and the shear viscosity, but increases the 
pressure; while the collapse is briefly accelerated. It is 
interesting to note that after the gaussian emission the 
distribution recovers staticity slowly, probably in a qua- 
sistatic regime. In this model pt > Pr, which means 
2V3aT] >0 {h^ 0.99). It is important to mention that 
in this model, a shear viscosity 77 > is only possible if 
we choose the negative root in p^ . Physically mean- 
ingful values of shear viscosity (77 > 0) are obtained for 



regions r/a « 0.6 — >■ 1. This means that the inner core 
is not viscous but anisotropic. The rest of the kinematic 
variables (O and a), shown in Figures 7-8, follow the 
evolution of the radius of the distribution, with Q (a) 
decreasing (increasing) faster as the radius decreases at 
a faster rate for 4 ;S t ;S 6. 
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FIG. 13: Evolution of the radii for the Tolman Vl-like model 
III. The initial conditions are A{0) = 8.0, -F(O) = 0.75, 0(0) = 
-0.1. 
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FIG. 14: Evolution of the local radial velocity uu for the Tol- 
man Vl-like model III The initial conditions are A{0) = 8.0, 
F(0) = 0.75, n{0) = -0.1, with h = 0.95. 



B. Schwarzschild— like model II: adiabatic 

We construct this model with the same effective vari- 
ables and metric functions as the aforestudied model I, 
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FIG. 15: Evolution of the shear viscosity 77 (multiphed by 
10'') for the Tolman Vl-like model III. The initial conditions 
are A{0) = 8, F{0) = 0.75, ^(0) = -0.1, with h = 0.95. 



distribution the inner shells collapse faster. The fffec- 
tive gravitation is therefore enhanced by the anisotropy 
induced by the viscosity. Inner regions have a greater 
shear viscosity in this model 10-10^ times the values 
found in model I). 



C. Tolman VI— like model III: non— adiabatic 

In this subsection we revise the model obtained from 
Tolman's solution VI [43 . Let us take 
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p 



gK{l - 9a{r/a)^^) 



3{K/I - 9a{r/a)^^^)hr'^ ' 
where K and / are defined as 



(43) 



(44) 



K ^8-3h + 4\/4 - 3/1, 



but now the radiation flux is zero everyhere; therefore 
this model is adiabatic. Obviously we do not need now 
an equation of state because all physical variables are 
determined algebraically from the field equations. How- 
ever some measure of tangential stress at the surface is 
required to evolve the system. We opt for a tangential 
pressure equal to the radial pressure just at the surface, 
Pt\a = Pr\a- The third surface equation in this case is 



hF - 6n*F + Shfl^F 



F - 3hn'^ - 4hn^ ~3n^ + i-h + en'^ 



(42) 



Observe that this expression explicitly depends on the 
anisotropic parameter h. In this case we integrate the 
system for the initial conditions ^(0) = 5, M(0) = 1 and 
r2 = —0.01. Figure 10 shows the radius of the distribu- 
tion for different values of h. Figures 11-12 display the 
radial velocity and the shear viscosity for h = 0.9. In 
this case anisotropy manifests clearly at the surface. As 
long as pt is greater than pr the collapse accelerates. The 
same occurs for 0.7 < h < 1.0, as we go deeper in the 



I = 8-3h- 4:^/4: - 3h. 

g and a are functions of which can be determined using 
(1361). Therefore 



9 = 



3(l-£) 
247r 



a = 



3h{l ~F)-KI3 



9[3/i(l - F) - 113] 
where 

p = {1 - F)n^ + 2E{i + n). 

Using (ini) and (HH) we obtain 

m = ma{r/a) 

and 



(45) 



(46) 



(47) 



V = InF H IH ln(r/a) H -1^= < / In ^ , ' . 

\^ a^^=^{K/I -9a) j] ' 

Evaluating equation (|23p at the surface we can obtain an equation for Q (too long to display here). 

Integrating the system of equations at the surface for the initial conditions A{0) — 8 and = —0.1, with Mr = 10~^, 
we obtain figures 13, 14 and 15. We obtain similar results as in model II but in a different fashion. The bigger the 
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difference between the tangential (viscous) pressure pt and the radial pressure Pr (as h decreases), the more violently 
the distribution explodes. It is striking that now all the spherical shells tend to reach the same instantaneous local 
radial velocity when the system goes to faster collapse with emission of energy across de boundary surface. At least 
locally, the "acceleration" of all the shells goes to zero at the same time; again the same instantaneous local radial 
velocity (negative) is reached before a final bouncing per shell from outer to inner. 



D. Tolman VI— like model IV: adiabatic 



We construct this model with the same effective variables and metric functions as in model III, but now the radiation 
flux is zero everyhere; therefore this model is adiabatic. Obviously we do not need now an equation of state because 
all physical variables are determined algebraically from the field equations. However some measure of tangential stress 
at the surface is required to evolve the system. We opt for a tangential pressure equal to the radial pressure just at 
the surface, pt\a = Pr|a, as in model II. The third surface equation in this case is 

tl = {-m'^Kh^I^F + fl'^FI^K^ + 9fl'^K^h^IF + GQ'^ F l'^ hy/ A - 3/i - 2A3FIh'^n'^ ~ I^K'^Ffl'^ 
- SlFh^I + MFIKh^^/A - ih + 2Amh'^Q^F - ISfl^ F l"^ K V A - 3h - ISQ"^ F I V A - 3h 

+ siFh'^K ~ mh'^i'^FKn^ + isFh^K^n^i + 9n^Kh^i^ - gn'^K^h'^i - ish^K^n^i + ish'^i^xn^ 

+ I^K^n"^ ~81h^K - I'^n-^K^ + 8lh'^I + 8lKh^n^ ~81Ih'^n^)/[l62{K ~ I)h'^A] (49) 



Integrating the system for the initial conditions A{0) = 8 
and = —0.1, figure 13 shows the radius of the distri- 
bution for different values of h. Figures 16-18 display 
the radial velocity and the shear viscosity for h = 0.9. 
After some numerical experimentation some non-trivial 
results arise, and we relax the condition pt > Pr- At the 
surface we do not find any novelty. The most violent ex- 
plosion occurs as » Pt- In this adiabatic but viscous 
(anisotropic) model all the shells bounce at the same time 
to irrupt from inner regions to outer regions with an ap- 
parently linear dependence with time. The outer shells 
of matter are ejected faster and earlier than the inner 
ones. This sort of behavior was reported several years 
ago studying in Bondi coordinates the collapse of radiat- 
ing distributions with an extreme transport mechanism 
as diffusion [46,]. However, the shear viscosity profiles in- 
dicate that i) bouncing is not allowed at all and ii) some 
inner regions are forbidden, otherwise the shear viscos- 
ity profiles become negative or/and infinite (see Figure 
18). This situation is general and independent of the 
anisotropy parameter h. 
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FIG. 16: Evolution of the radii for the Tolman Vl-like model 
IV. The initial conditions are A{0) = 8.0, F{0) = 0.75, 0(0) = 
-0.1. 



V. CONCLUSIONS 

We consider a selfgravitating spherical distribution of 
matter containing a dissipative fluid. The use of the 
PQSA with non-comoving coordinates allow us to study 
viscous fluid spheres in the streaming out limit as they 
just depart from equilibrium. From this point of view, 
the PQSA can also be seen as a nonlinear perturbative 
method to test the stability of solutions in equilibrium. 

For the non-adiabatic Schwarzschild model the distri- 
bution evolves to a flnal state with a non-viscous and 



anisotropic inner core. Surprisingly, in this model the 
evolution of the local radial velocity at the surface is 
the same in the isotropic (pt — Pr) fortuitous 
coincidence due to the chosen equation and state and 
the incompressibility of the fluid. For the adiabatic 
Schwarzschild model the final core is up to 10^ times 
more viscous, and the anisotropy appears explicitly in all 
the evolution equations. The higher viscosity of the core 
increases the effective gravity and the collapse is faster, 
as long as pt > Pr- 

Both of the Tolman VI models lead to a distribu- 
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FIG. 17: Evolution of the local radial velocity u for the Tol- 
man Vl-like model IV. The initial conditions are ^1(0) = 8.0, 
F{0) = 0.75, n(0) = -0.1, with h = 0.95. 
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tion which initially collapses and then bounces and 
expands indefinitely. The Tolman VI non-adiabatic 
model shares some of the characteristics of the adiabatic 
Schwarzschild. Before the final bouncing, as pt > Pr the 
collapse is accelerated. For the non-adiabatic case some 
regions of the parameter space are forbidden, since the 
shear viscosity profiles become unphysical. In this case 
the bouncing is not allowed and the distribution collapses 
indefinitely. 

A forthcoming paper considers the dissipation by heat 
flow, in order to isolate effects similar to the ones stud- 
ied in the present investigation, but with different mech- 
anisms. Also, a work in progress considers heat flow and 
anisotropy induced by electric charge, pointingto the 
most realistic numeric modeling in this area j47l |. Al- 
though they are not entirely new, the results presented 
here constitute a first cut to more general situations us- 
ing the PQSA, including dissipation, anisotropy, electric 
charge, heat flow, viscosity, radiation flux, superficial ten- 
sion, temperature profiles and study their influence on 
the gravitational collapse. This investigation is an essen- 
tial part of a long-term project which tries to incorporate 
the Miiller-Israel-Stewart theory for dissipation and de- 
viations from spherical symmetry, specially when con- 
sidering electrically charged distributions. Besides being 
interesting in their own right, we believe that spherically 
symmetric fluid models are useful as a test bed for more 
general solvers in numerical relativity [39^, '40'] . A general 
3D code must be able to reproduce situations closer to 
equilibrium. 
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